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ABSTRACT 

Marked correlations are reported between the ionization of the X-ray line 
spectra of normal OB stars, as observed by the Chandra X-Ray Observatory, and 
their optical spectral types. These correlations include the progressive weakening 
of the higher ionization relative to the lower ionization X-ray lines with advancing 
spectral type, and the similarly decreasing intensity ratios of the H-like to He- 
like lines of the a ions. These relationships were not predicted by models, nor 
have they been clearly evident in astrophysical studies of a few objects; rather, 
they have emerged from morphological analysis of an adequate (albeit still small) 
sample, from which known peculiar objects such as magnetic stars and very rapid 
rotators have been isolated to reveal the normal trends. This process is analogous 
to that which first demonstrated the strong relationships between the UV wind 
profiles and the optical spectral types of normal OB stars, which likely bear a 
physical as well as a historical connection to the present X-ray results. Since the 
optical spectral types are calibrated in terms of fundamental stellar parameters, 
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it follows that the winds and X-ray spectra are determined by the latter. These 
observations provide strong guidance for further astrophysical modeling of these 
phenomena. 

Subject headings: stars: early-type — stars: fundamental parameters — stars: 
winds — X-rays: stars 

1. Introduction 

When a new observational domain, comprising a diverse phenomenology not yet fully 
explained or predicted by physical models, is opened to investigation, the science of morphol- 
ogy may provide a useful initial formulation. For instance, morphology can discern empirical 
trends, separate the normal (majority) behavior from peculiar exceptions, reveal correlations 
with other kinds of information, and suggest or eliminate hypotheses for subsequent physical 
analysis, all with a minimum of prior assumptions. A frequent occurrence of this situation 
in astronomy is access to a new wavelength region or data of higher information content. A 
specific example from the relatively recent past is the extensive database of high-resolution 
ultraviolet stellar spectra generated by the International Ultraviolet Explorer (IUE) satellite. 
During the early 1980s, there was a debate about the relationship if any of OB stellar winds 
to the fundamental stellar parameters. A majority of specialists believed that it was weak 
or nonexistent, because of the large scatter in plots of mass- loss rate versus effective temper- 
ature or luminosity; however, it is now known that the derivations of those quantities were 
inadequate at that time. The debate was settled by a purely morphological demonstration 
of the detailed correlations between the wind profiles in the IUE data and the optical spec- 
tral types, which had in turn been calibrated in terms of fundamental parameters (Walborn 
et al. 1985, 1995, 2009; Rountree & Sonneborn 1991, 1993; Penny et al. 1996). The database 
provided by the Far Ultraviolet Spectroscopic Explorer (FUSE) supported the extension of 
those correlations to the rich spectral region below Lyman-a (Walborn et al. 2002, Pellerin 
et al. 2002). Of course, morphology does not explain anything, but it may pave the way 
toward the ultimate objective of physical understanding. 

The Chandra (and XMM-Newton) X-ray satellites have provided the first stellar spec- 
troscopic data of adequate quality to investigate the systematics of that regime. The samples 
remain small to date because the required high-resolution data are time-consuming to ac- 
quire with those systems. Nevertheless, they are sufficient for an initial morphological survey, 
which is presented here based on the Chandra data, with results that may be surprising. 
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2. Observations 

All of the observations were obtained with the Chandra Advanced CCD Imaging Spec- 
trometer (ACIS) in its High Energy Transmission Grating Spectrometer (HETGS) config- 
uration, which in turn simultaneously provides Medium Energy Grating (MEG) and High 
Energy Grating (HEG) data. The resolving powers are 660 at 15 A and 1000 at 12.4 A. We 
have combined the two datasets whenever possible, i.e., whenever the HEG S/N was ade- 
quate. The data discussed here were obtained from the TGCat (Huenemoerder et al. 2009) 
web-based catalog of Chandra grating data. TGCat processing starts with the Level 1 event 
data downloaded from the Chandra archive. Updated calibrations and responses are applied 
to the data and the zeroth order position is determined to high accuracy. The resulting 
TGCat Level 2 extracted spectral data were used to create the plots shown here. 

Most of the data discussed here were archival, but three of the four targets observed 
in our Chandra program (PI Waldron, ID 6200204) to fill gaps in the extant coverage of 
the HR Diagram are included. They are HD 93250, 9 Sgr, and 15 Mon; our fourth target, 
HD 93129AB, presents a complex, composite spectrum that will be discussed separately 
(Nichols et al., in preparation.) Stellar and observational parameters are listed in Table 1, 
in the order in which the objects appear in Figures 1-6, together with extensive references 
to previous analyses of these stars. Further explanation of several columns in the table is 
given below. 

Other relevant data exist in the XMM-Newton archive, which we plan to explore. Those 
observations provide better coverage of the CNO lines at the longer wavelengths, and they 
include several stars not observed by Chandra. 

3. Results 

The 5-25 A X-ray spectra of 14 normal OB stars (or binary systems) are shown in 
order of advancing optical spectral types in Figures 1 and 2, which contain primarily main- 
sequence stars and supergiants, respectively, albeit with a few giants included to augment 
the coverage in each case. In these plots, each X-ray spectrogram has been scaled to 1.0 at 
the peak emission feature. Thus, absolute intensities cannot be directly compared between 
different objects, although line ratios can be. Because of the large range of line intensities in 
the calibrated data, a complementary presentation of the same data, but with three different 
scalings at wavelengths below 11 A, from 11 to 18 A, and above 18 A is given in Figures 3 
and 4. The three normalization factors for each object are listed in Table 1; their units are 
photons s" 1 cm -2 A -1 . The largest of the three values (usually the third, except for three 
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of the earliest types where it is the second) applies to Figures 1 and 2. As indicated by the 
spectral types, several objects are known close binaries. Parameters of the strong spectral 
lines are listed in Table 2. 

Several correlations between the X-ray line ionization and the optical spectral types 
are immediately apparent in these figures. First, in Figures 1 and 2, the higher ionization, 
shorter wavelength lines weaken markedly relative to the lower ionization, longer wavelength 
lines in a given spectrum along the sequences, although the oxygen lines are usually the 
strongest. Second, in Figures 3 and 4, the behaviors of individual lines and line ratios, 
particularly those of the close pairs of lines from the H-like and He-like a ions, display 
strong trends along the sequences. Specifically, Sxv essentially disappears after 03.5 on 
the main sequence, and after 04 in the supergiants. The Sixiv/Sixm and Mgxn/Mgxi 
ratios diminish drastically between 03.5 and 04 in both sequences. The Nex/Neix ratio 
is the best diagnostic, because it is the most sensitive throughout this spectral-type range. 
Unfortunately, the Nex feature is blended with an Fexvil line that strengthens toward the 
later types (as do the other lines from the latter ion), interfering with the visual trend; 
this blend is elucidated at a larger scale in Figure 5. Finally, the Ovm/Ovil trend is 
relatively weak throughout this range, although the ratio reverses between the extreme cases 
of HD 150136 and f3 Cru (the spectrogram of HD 93250 has lower S/N and is unreliable at 
the longer wavelengths). The large N/O ratios in ( Pup and £ Per are most likely abundance 
effects in processed material, as previously derived in the former spectrum by Kahn et al. 
(2001) in agreement with prior photospheric/wind analyses; see also Oskinova et al. (2006). 
Note that the differential behaviors of the ionization ratios are themselves an ionization 
effect, with the intermediate Ne lines providing the best response to the relevant parameter 
range. Quantitative measurements (and corresponding uncertainties) of these line ratios, as 
well as the implied temperatures, are given by Waldron & Cassinelli (2007). We present 
measurements of the Ne ratios below. 

It is noteworthy that the binary nature of several objects, which may have colliding 
winds, does not interfere with these observed trends, although it may contribute to some 
of the scatter in them. Neither does the light to moderate range of interstellar extinctions 
among them (Table 1), for which no correction has been made here. The seven stars with 
E(B-V) between 0.3 and 0.5 are expected to have emission features longward of 18 A extincted 
by factors up to 2-3. One might be concerned that the earlier type normal stars tend to 
have the higher interstellar extinctions, but ( Pup provides the counterexample with an 
early type and low extinction that fits the X-ray line sequence well. Of course, wind or 
interstellar extinction always affects the longer wavelengths more, so it cannot cause the 
relatively weaker high ionization features at the later types. Neither would it affect the 
ratios of the close pairs of H-like to He-like lines significantly. 
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The present sample is not adequate to investigate luminosity effects in detail, but there 
are a few suggestive indications of an ionization correlation with that dimension as well. 
As noted above, the Sxv line persists to a slightly later type in the supergiants. More 
systematically, it can be seen that the Si xiv and Mg XII lines remain stronger at later types 
in the supergiants than on the main sequence (see also Waldron & Cassinelli 2007). Indeed, 
i Ori rather breaks these trends in Figure 3, which may well be due to its giant nature. 
More extensive coverage of intermediate luminosity classes with X-ray data of this quality 
is required to determine whether these effects correlate in detail with that dimension. 

Figure 6 displays the X-ray spectra of two magnetic and two very rapidly rotating OB 
stars. It is easily seen that they have much higher ionizations than those of normal stars with 
the same spectral types. They undoubtedly correspond to energetic circumstellar activity 
related to the fields or rotation (see also Wojdowski & Schulz 2005, Zhekov & Palla 2007, 
Westbrook et al. 2008)0 This result emphasizes once again that it is essential to eliminate 
pathological cases from the samples, if one intends to investigate the normal phenomena. 

To summarize and somewhat quantify our results, we have measured the optimal Ne 
ratio, by means of Gaussian fits to the lines. Since the Nex blend with the Fexvil 12.124 A 
line cannot be resolved, we have corrected for it approximately by measuring the isolated 
Fexvil 15.014 A line and subtracting it, scaled by the theoretical intensity ratio of 0.13, 
from the blend. The resultant Ne ratios are listed in Table 1 and plotted versus spectral 
type in Figure 7, with the luminosity classes and the peculiar objects distinguished by 
different symbols; the trends described above are well shown. We have also derived X-ray 
temperatures from the Ne ratios by procedures similar to those of Waldron & Cassinelli 



(2007)13 they are likewise listed in Table 1 and plotted in Figure 7 against the stellar effective 
temperatures from the recalibration of the spectral types by Martins, Schaerer, & Hillier 



*As shown in Table 1, £ Oph and 7 Cas have substantially higher vsini than the other stars in this sample. 
However, they are likely not typical rapid rotators, and it should not be inferred that all objects in that 
category would have anomalous X-ray spectra. £ Oph displays transient Oe/Ha emission events (Niemela & 
Mendez 1974, Ebbets 1981), while 7 Cas is subject to extreme shell episodes (Cowley & Marlborough 1968) 
and may also be magnetic (Smith et al. 2004). 

2 It should be noted that the present Nex/Neix ratios are systematically smaller than those of Waldron 
& Cassinelli (2007), for two main reasons: they chose not to remove the Fexvil line from the Nex blend, as 
was done here; and they estimated the Ne IX r-component in isolation from the f, i and other blended species, 
whereas we measured the entire complex. However, the respective theoretical determinations approximate 
the different contents of the measurements, so that the temperatures derived from the two studies agree within 
±16% in all cases except that of HD 206267 (35%). In this spectrum, the Fexvil line at 15.014 A, which 
we use to estimate the strength of the Fexvil 12.124 Aline blended with Nex, appears to be anomalously 
weak, leading to a possible undercorrection and hence Nex overestimate here (Waldron, in preparation). 
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(2005). The two plots are very similar, as expected from the good correlations between the 
respective variables. No systematic dependence of the Ne X-ray ionization on the luminosity 
classes is evident in this sample, nor does a plot against the gravities show any meaningful 
trend (e.g., the normal class V stars span the full X-ray ionization range but have constant 
gravity) . 

4. Discussion 

The systematic characteristics of OB X-ray spectra just described have not emerged 
clearly from previous studies. The likely reasons are inadequate samples contaminated by 
peculiar objects, and the uncertainties of current physical models for the phenomena. As 
recounted in the Introduction, the situation regarding the systematics of OB stellar winds 
during the early 1980s was remarkably similar, and it was substantially illuminated by purely 
morphological investigations. In fact, the mass- loss rates and physical details of OB winds 
remain uncertain today, with the roles of clumping, rotation, and possibly weak magnetic 
fields being actively investigated. Nevertheless, the purely empirical descriptions of the 
phenomena stand as guides and objectives for the modeling. 

The facts that the UV wind and X-ray spectral morphologies both correlate with the 
optical spectral types is likely no coincidence. Indeed, the X-rays most likely originate 
in the winds (Lucy & White 1980, Owocki et al. 1988, Feldmeier et al. 1997), and the 
observed correlations indicate that the detailed wind structures must be well determined 
by the fundamental stellar parameters, likely via the intense radiation fields of the OB 
stars. Waldron & Cassinelli (2007) measured and analyzed several parameters in the X-ray 
spectra of nearly the same sample discussed here. The physical details are very complex, 
e.g., the lines considered here must form at different depths in the winds and are thus 
subject to differing amounts of internal extinction. X-ray emission from shocks related to 
clumps distributed throughout the winds is discussed by Feldmeier et al. 2003, Oskinova 
et al. (2004, 2006), Zhekov k Palla (2007), Cassinelli et al. (2008), and Waldron & Cassinelli 
(2009). There may also be global asymmetric wind structures even in the normal stars, e.g., 
two-component equatorial vs. polar winds related to rotation and/or weak magnetic fields 
(Mullan & Waldron 2006). Nevertheless, the global systematics evident in the emergent 
spectra strongly suggest that all of these physical effects must be correlated and determined 
by the stellar and wind parameters. This result in turn provides strong guidance for further 
physical modeling of the winds and X-ray generation. 

This conclusion may appear surprising, yet it was foreshadowed by the strong correlation 
between the bolometric and X-ray luminosities of the OB stars that has been known for many 
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years (Long & White 1980, Sana et al. 2006). Moreover, several of the objects with peculiar 
X-ray spectra, such as those in Figure 6, deviated from that correlation (e.g., Cassinelli et al. 
1994, Cohen et al. 1997). Of course, further physical modeling that uniquely reproduces the 
morphological correlations is required in order to understand them. 
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WLW also acknowledges partial support from SAO grants AR8-9003A and GO8-9021B. 
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Figures 2 and 4: 
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3.46 


0.69 


Waldron 


17, 18, 


19, 20 


e Ori 


BO la 


30.0 


91 


0.08 


4.15-4/3.28-3/8.25-3 


0.22 


3.46 


1.06 


Waldron 


21 





Figure 6: 



e 1 Ori c 


07 Vp 


36.9 


130 


0.32 


3.68-3 


0.87 


5.04 


4.59 


Canizares, Gagnc 


22. 


23 


T SCO 


B0. 2 V 


30.0 


<5 


0.04 


1.23-2 


0.42 


4.08 


0.83 


Cohen 


24, 


25 


C Oph 


09.5 Vnn 


31.9 


385 


0.32 


4.00-3 


0.43 


4.11 


0.97 


Waldron 


26, 


27 


7 Cas 


B0. 5 IVc var 


30.0 


380 


0.08 


5.72-3 


1.64 


6.23 


0.60 


Smith 


28 





"Martins et al. 2005; 30 kK adopted for B0-B0.5. 

b Howarth & Prinja 1989, Howarth ct al. 1997, Harmanec 2002 (7 Cas), SIMBAD (f) Cru). 

c Pcak emission fluxes in photons s" 1 cm -2 A -1 by which the spectrograms have been divided to normalize them in the plots. The negative powers of ten follow the minus 
signs. The first value pertains to wavelengths below 11 A, the second for 11—18 A, and the third above 18 A in Figures 3 and 4; whereas, only the largest (or sole) value applies 
to the entire wavelength ranges in Figures 1, 2, (and 6). 

References. — References to Table 1: 1, Evans ct al. 2004; 2, Rauw ct al. 2002; 3, Burkholder et al. 1997; 4, Gies et al. 1997; 5, Gies et al. 1996; 6, Skinner et al. 2008; 7, 
Sanz-Forcada et al. 2004; 8, Cohen et al. 2008; 9, Skinner ct al. 2005; 10, Niemcla & Gamen 2005; 11, Lcutencgger ct al. 2007; 12, Kramer et al. 2003; 13, Cassinclli ct al. 2001; 
14, Kahn et al. 2001; 15, Harvin ct al. 2002; 16, Miller ct al. 2002; 17, Raassen et al. 2008; 18, Pollock 2007; 19, Cohen et al. 2006; 20, Waldron & Cassinelli 2001; 21, Waldron 
2005a; 22, Gagnc et al. 2005; 23, Schulz et al. 2003; 24, Cohen et al. 2003; 25, Mewe ct al. 2003; 26, Li et al. 2008; 27, Waldron 2005b; 28, Smith et al. 2004. 

Note. — The references in the table are to studies of individual stars, in inverse chronological order for a given star. Analyses of many or most stars in the table have also 
been presented by Waldron & Cassinclli 2007, Zhekov & Palla 2007, Wojdowski & Schulz 2005, and other references cited in the text. 



Table 2. Rest Wavelengths of X-Ray Lines [A] 
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Fig. 1. — X-ray spectral-type sequence of normal dwarfs and giants. The ordinate units 
are normalized to 1.0 at the peak emission feature in each spectrum, and the corresponding 
factors are given in Table 1. Thus, this scale may not be intercompared among different 
spectrograms. The wavelengths of the identified spectral lines are given in Table 2. 
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Fig. 2. — X-ray spectral-type sequence of normal supergiants and giants. The ordinate units 
and the wavelengths of the identified spectral lines are as in Fig. I. 
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Fig. 3. — The same X-ray spectral-type sequence of normal dwarfs and giants as in Fig. 1, 
except that the ordinate units are normalized to 1.0 by different factors in each of the three 
panels, at the peak emission feature in each. The corresponding factors are given in Table 1. 
The wavelengths of the identified spectral lines are given in Table 2. 
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Fig. 4. — The same X-ray spectral-type sequence of normal supergiants and giants as in 
Fig. 2, except with different normalization factors in each panel, as in Fig. 3. The wavelengths 
of the identified spectral lines are given in Table 2. 
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Fig. 5. — An expanded wavelength scale for the Nex feature, blended with an Fexvn line, 
in all the spectra of Figures 1-4. While the two features cannot be resolved, the increasing 
relative contribution of the Fe line at a slightly shorter wavelength (Table 2) in the later 
spectral types can be discerned. The Ne feature is very broad in the two earliest dwarfs, and 
it is blueshifted in the two earliest supergiants. 
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Fig. 6. — X-ray spectrograms of the magnetic stars 6 1 Ori C (Donati et al. 2002, Gagne 
et al. 2005, Wojdowski & Schulz 2005, Smith & Fullerton 2005, Stahl et al. 2008) and r Sco 
(Cohen et al. 2003, Wojdowski & Schulz 2005, Donati et al. 2006); and the rapid rotators 
C Oph (Waldron 2005b, Li et al. 2008) and 7 Cas (Smith et al. 2004). The H/He-like line 
ratios are much larger than in normal stars of the same spectral types. The ordinate units 
and the wavelengths of the identified spectral lines are as in Fig. 1. 
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Fig. 7. — Left: Plot of the Nex/Neix ratios as a function of spectral type. The symbols are 
as follows: filled squares, luminosity class V; filled triangles, classes III and II; filled circles, 
class I; and open circles, peculiar objects. The symbols for i Ori and r CMa have been 
offset by ±0.018 vertically, to avoid complete superposition. Right: Analogous plot of the 
X-ray temperatures derived from the Ne ratios vs. the effective temperatures corresponding 
to the spectral types (Table 1). t Ori and r CMa have been offset by ±0.025 vertically. In 
each plot, the normal stars show a declining X-ray ionization trend with advancing type or 
effective temperature, while the peculiar objects have higher ionization than normal stars of 
the same types or temperatures. 



